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Abstract. The stochastic dynamics of c and b quarks in the fireball created in 

nucleus-nucleus collisions at RHIC and LHC is studied employing a relativistic 

Langevin equation, based on a picture of multiple uncorrelated random collisions 

^. ■ with the medium. Heavy-quark transport coefficients are evaluated within a pQCD 

VO , approach, with a proper HTL resummation of medium effects for soft scatterings. The 

Langevin equation is embedded in a multi-step setup developed to study heavy-flavor 
f -) . observables in pp and AA collisions, starting from a NLO pQCD calculation of 

f " initial heavy-quark yields, complemented in the nuclear case by shadowing corrections, 

f~^ ■ fey-broadening and nuclear geometry effects. Then, only for AA collisions, the Langevin 

equation is solved numerically in a background medium described by relativistic 
hydrodynamics. Finally, the propagated heavy quarks are made hadronize and decay 
into electrons. Results for the nuclear modification factor Raa of heavy-flavor hadrons 
and electrons from their semi-leptonic decays are provided, both for RHIC and LHC 



X 



$H ' beam energies. 

Heavy-flavor electron spectra, measured in Au-Au collisions at y/s NN = 200 GeV by 
PHENIX |TJ |2] and STAR |3] experiments at RHIC, have displayed a large suppression 
with respect to pp collisions, comparable in amount to the one observed for charged 
hadrons. Models considering medium-induced gluon radiation [U [5] as the dominant 
energy loss mechanism for heavy quarks propagating in QGP come up against difficulties 
in reproducing such results. 

These findings gave a great boost to calculations taking in consideration the role 
of collisional energy loss [6j[7]. In some models the heavy-quark propagation in QGP is 
described through a Langevin stochastic equation [8j El [TO] , assuming that heavy-quark 
spectrum modifications arise from the cumulated effect of many uncorrelated random 
collisions with the medium. 



In our approach we use a relativistic Langevin equation [IT], describing the time 
evolution of the heavy-quark momentum: 

^ = -Mp)p t +e(t) (i) 

that involves a deterministic friction term and a stochastic noise term £*(£), completely 
determined by its two-point temporal correlator: 

(em j (t')) = b ij (p)6(t - £'), with b ij (p) = K L (pW + K T (p)(S ij - fpi) (2) 

The latter involves the transport coefficients kt{p) = \ ^ and kl{p) = J^+ , which are 
evaluated according to the procedure presented in [12] . We introduce an intermediate 
cutoff \t\* ~m|) (t = (P'—P) 2 ) to separate hard and soft scatterings. The contribution of 
hard collisions (|t| > \t\*) is evaluated through a pQCD calculation of the processes 
Q(P)qi/i —> Q{P')qi/i and Q(P)g — > Q(P')g. On the other hand for soft collisions 
(|t| < |t|*) a resummation of medium effects is provided by the Hard Thermal Loop 
approximation, with a s (/x) evaluated at a scale fi oc T. The final result [12j is given by 
kt/l(p) = ^t/lG ) + k t/l(p)- According to the scale \i at which a s (fi) is evaluated to 
calculate K^jl-, we devised two different sets of calculations, referred to in the following 
as HTL1 (for \i oc T) and HTL2 (for fi = \t\). 

The Langevin simulation tool is embedded in a full setup to calculate heavy-flavor 
observables in pp and AA collisions, divided into the following independent steps [12] : 

(i) A sample of c and b quarks is generated using POWHEG [13], a code which 
implements pQCD at NLO accuracy, with CTEQ6M PDFs as input. For AA 
collisions, EPS09 nuclear corrections to PDFs are employed [H]; then, heavy 
quarks are distributed in the transverse plane according to the nuclear overlap 
function Tab(x, y) =T^(x+6/2, i/)Tb(x— 6/2, y) corresponding to the selected impact 
parameter b; a p^-broadening correction to heavy-quark momenta is also included. 

(ii) At a given proper-time r an iterative procedure is started, only for AA collisions, 
to follow the stochastic evolution of the heavy quarks in the plasma until 
hadronization: the Langevin transport coefficients are evaluated at each step 
according to the local 4- velocity and temperature T(x) of the expanding background 
medium, as provided by hydrodynamic codes: both ideal and viscous fluid scenarios 
were tested [T5] [T6l [T7] . 

(iii) Heavy quarks are made fragment into hadron species, sampled according to c and 
b branching fractions taken from Refs. [THJ [19] , while their momenta are sampled 
from a Peterson fragmentation function [20], with e = 0.04 and 0.005 for c and 
b respectively. Finally, each heavy-quark hadron is forced to decay into electrons 
with PYTHIA [21], using updated tables of branching ratios based on Ref. [19J. 

The effects of the Langevin evolution of c and b quarks in AA collisions resulting 
from our calculations are studied through the nuclear modification factor Raa{.Pt) = 
{dN/dp T ) AA /(N) con {dN/dp T ) pp and the elliptic flow coefficient v 2 {p T ) = (cos(2(/>)) PT of 
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Figure 1. (left) Raa(pt) of heavy-flavor electrons in minimum-bias Au-Au collisions 
at RHIC (0 — 92% of the inelastic cross section) ; (right) Raa obtained by integrating 
the electron yields over the indicated momentum ranges, as a function of N part . In both 
panels our predictions (with HTL1 calculation) are compared to PHENIX data [2]. 



final-state heavy-quark hadrons or decay-electrons. Here we will display only results 
obtained using viscous hydrodynamics for some representative values of the input 
parameters (r and the QCD scale /x), among those fully explored in [12] . 

In Fig. [I] our findings for the Raa of heavy-flavor electrons in Au-Au collisions at 
RHIC (y/s NN = 200 GeV), obtained with the calculation HTL1 by assuming viscous 
hydrodynamics, r = 1 fm/c and /i = 3nT/2, are compared to PHENIX [2] data. In 
the left panel we observe a general agreement of Raa{pt) with the PHENIX results on 
a minimum-bias data sample (0 — 92% of the inelastic cross section) for pr ~ 3 GeV/c. 
In the right panel, displaying the Raa obtained by integrating the electron yields above 
a given pt and plotted versus N part , the centrality dependence of the PHENIX data is 
shown to be nicely reproduced, except for the most-peripheral centrality bin. We do not 
show here the results obtained for the heavy-flavor electron t>2 (pt) [12], that appear to 
underestimate the PHENIX minimum-bias data. However, a more detailed treatment 
of hadronization, including also the coalescence mechanism, should enhance both v-i and 
Raa at pt ~ 3 GeV/c. 

In Fig. [2] we show our predictions for the Raa{pt) of heavy-flavor electrons (left 
panel) and of D, B mesons (right panel) in Pb-Pb collisions at LHC (y/s NN = 2.76 TeV), 
calculated by selecting the — 10% most-central events, under the hypothesis of viscous 
hydrodynamics and for two different choices (HTL1 or HTL2) of the \x scale in the 
calculation of K^T^. General features of the Raa of heavy-flavor electrons appear similar 
to those observed at RHIC at the same centrality |12j . with a stronger suppression of 
both charm and bottom contributions. However, results obtained with the calculation 
HTL2 display a weaker quenching for p T > 3 GeV/c. As regards D and B suppression, 
slighter differences between HTL1 and HTL2 are observed, and at higher pt- 

The main goal of our study was to deliver a weak-coupling calculation to be used 
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Figure 2. The Raa of heavy-flavor electrons (left) and D, B mesons (right) in Pb-Pb 
collisions at LHC (y/s NN = 2-76 TeV) for the — 10% most-central events obtained 
with the HTL1 (solid line) and HTL2 (dashed line) calculations. 



as a benchmark for advanced studies or less conventional scenarios. The capability to 
accommodate the electron spectra observed at RHIC for p T ^ 3 GeV/c strengthens the 
hypothesis that heavy-quark collisional energy loss must be taken into proper account. 
Moreover, we will soon be able to test our predictions against first data with Pb-Pb 
collisions delivered by LHC experiments [22j [231 Wi 
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